Microwave power generated by a spin-torque oscillator in the presence of noise 
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An expression for the microwave power generated by a spin-torque oscillator in the presence of 
thermal noise is derived. This expression, when used in a subcritical (I < / t h) regime, demonstrates 
that generated power is determined by the noise level and provides a simple way to experimentally 
determine the threshold current 7 t h for microwave generation. The same expression gives a good 
quantitative description of the experimentally measured dependence of the generated power on the 
bias current p(I) in a moderate super-critical (ith < I < 1.57th) range of bias current variation. 



It has been demonstrated both theoretically [J and 
experimentally [1,0, Q that the spin-torque effect in mag- 
netic nanostructures creates an effective negative mag- 
netic damping, which for sufficiently large bias currents 
(i > 2th) can compensate positive damping. This leads 
to the excitation of one of the spin wave eigenmodes of 
the structure, thus providing a way to create fully metal- 
lic nano-sized microwave spin-torque oscillators (STO). 

In contrast with almost all other conventional mi- 
crowave oscillators, STOs are very small in size and, 
therefore, the influence of thermal fluctuations (thermal 
noise) on their operation is of a critical importance. This 
problem bears some similarity with the case of microwave 
oscillators based of Josephson junctions @ ■ An STO in 
the sub-critical regime (i < ith) in the presence of ther- 
mal noise, in particular, functions as a regenerative filter. 
With the increase of bias current in this regime (which 
corresponds to a decrease of net magnetic damping), the 
oscillator filters effectively its lowest eigenfrequency from 
noise, thus providing a measurable output microwave sig- 
nal at this frequency much before the threshold of self- 
sustained microwave generation (I = ith) is attained. As 
such, an experimental determination of the STO thresh- 
old current in the presence of thermal noise presents a 
challenge, which can only be met if the influence of ther- 
mal noise on the STO operation is accounted for in the- 
ory. 

In this Letter we present a theory of microwave gener- 
ation in STOs in which the influence of thermal noise is 
explicitly taken into account, in contrast with previously 
developed theories 0, [|| ■ The developed spin- wave the- 
ory of microwave generation is valid for a wide range of 
bias currents: from zero current at which the system is in 
the state of thermal equilibrium, to a significantly above- 
threshold (super-critical) currents at which the influence 
of thermal noise can be neglected. The theory gives good 
quantitative description of recent experimental results on 
microwave power generated by a current-driven magnetic 
nanopillar [sj, and also allows for a simple and accurate 
procedure for determining the threshold current of self- 
sustained microwave generation in STOs from the mea- 
sured dependence of the output microwave power on the 
bias current in the subcritical regime. 



It has been demonstrated in Refs.0, Hthat the magne- 
tization dynamics in STOs can be adequately described 
by the dimensionless complex amplitude c(t) of a sin- 
gle spin wave mode excited by the bias direct current. 
The phase <f> = arg(c) of the complex amplitude c(t) is 
equal to the azimuthal angle of the magnetization pre- 
cession, while the oscillation power p = \c\ 2 determines 
the polar angle of the precession, 9 = arccos(M z /Mo) = 
arccos(l — 2p). Here Mq is the length of the magnetiza- 
tion vector in the "free" magnetic layer, and M z is the 
projection of this vector on the direction of stationary 
equilibrium magnetization z (see Ref. for details) . 

The equation of motion for the amplitude c(t) in the 
general case, whereby the thermal noise is taken into ac- 
count, can be written as 



dc 
di 



+ iuj(p)c + r + (p)c-T_(p)c = f(t), 



(1) 



where u)(p) is the frequency of the excited spin wave 
mode, r + (p) is the natural positive damping in magnetic 
system, T_(p) is the effective negative damping intro- 
duced into the system by the spin-polarized current, and 
fit) is a random white Gaussian process that describes 
the influence of the thermal fluctuations (noise) . The cor- 
relation function of this random noise f(t) can be written 
as 



(f*(t)f(t')) = 2DS(t-t') , 



(2) 



where D is the diffusion coefficient that characterizes the 
noise amplitude. 

In general, the frequency Lo(p) and the effective damp- 
ing rates T±(p) arc complicated nonlinear functions of 
the spin wave power p. For moderate levels of excita- 
tions of the spin system, however, one can use a Taylor 
series expansion of these functions and write approxi- 
mately (see, e.g., Refs. ITol andflll): 



r + (p)«r (i 



Lo(p) w ujQ + Np, (3a) 
Qp), r_(p)w(rJ(l-p).(3b) 



Here loq is the linear oscillation frequency and N is 
the nonlinear frequency shift coefficient (see Eq. (38) in 
Ref. 0), Tq is the linear damping of the excited spin wave 



2 



mode in the passive regime, Q is a phenomenological coef- 
ficient characterizing the nonlinearity of positive damp- 
ing (see Ref. [Ill for details), a is the spin-polarization 
efficiency defined in Eq. (2) of Ref. [l(J and / is the bias 
charge current. 

We would like to stress that Eq. can adequately 
describe an auto-oscillator of any nature under the influ- 
ence of white noise fit), provided that this oscillator has 
a nonlinear frequency uj(p), nonlinear natural positive 
damping r + (p), and nonlinear negative damping r_(p), 
where p = |c| 2 is the oscillation power. 

It should also be noted that the complex spin wave am- 
plitude c is a canonical complex variable of the magnetic 
system. In particular, this means that the phase volume 
dfl of some element of the system's phase space is simply 
proportional to the area of this element on the complex 
c-plane, i.e., in power-phase coordinates, d£l = dpd(f>. As 
such, the probability distribution function V cq (p,4>) of 
the system at thermal equilibrium (i.e. for r_(p) = 0) 
should have the form of the pure Boltzmann exponential 
distribution, 



V eq (p, (f>) oc exp 



k B Tj 



(4) 



without any pre-exponential factors. In Eq. Q k B is the 
Boltzmann constant, T is the absolute temperature, and 
£ (p) is the energy of the system, 



£{p) = A / Lu(p')dp' 



A j uj p+ T^Np 2 



(5) 



Here the constant A depends on the normalization of the 
amplitude c and, in the case of an STO, is given by 



A = Fe ff Mo/7 , 



(0) 



where V c g is the effective volume of the magnetic material 
of the "free" layer involved in the auto-oscillation, and 7 
is the gyromagnetic ratio. 

To describe correctly the stochastic dynamics of a non- 
linear oscillator with arbitrary dependences of the fre- 
quency lo = Lu(p) and natural damping T + = r + (p) on 
the oscillation power p, one has to assume that the dif- 
fusion coefficient D also depends on p as 



D(p) = r + (p) v (p) = r + ( P ) 



Xu(p) 



(7) 



where 77 (p) is the effective noise power in the nonlinear 
regime. As it will be shown below, this form of the diffu- 
sion coefficient provides correct equilibrium distribution 
Eq. Q for arbitrary power dependences uj(p) and T + (p). 

The stochastic Langevin Eq. |T]) produces the following 
deterministic Fokker-Planck equation for the probability 
distribution function V(t,p, (f>) of the system, 



dV 
~dt 



d_ 

dp 



[2 P (r + -r_)p] 



dV 



(8) 



( 2p D—\ + —(— — 

dp \ dp J d4> \2p dej) 



While the general non-stationary solution of the 
Fokker-Planck Eq. is rather complicated, it is possi- 
ble to find a simple analytic expression for the stationary 
probability distribution function 'Po(p) that is indepen- 
dent of the oscillator phase 0, since in the stationary state 
all the phases of the oscillator are equally possible. This 
expression has the form, 



Va(p) = A^exp 



A 



W) (1- 



r-(pQ 



dp' 



where A/"o is the normalization constant determined from 
the normalization condition J °° Vo(p)dp = 1. 

From Eq. ([9]) one can see by inspection that our ap- 
proach gives the standard Boltzmann distribution Eq. 
at thermal equilibrium (r_(p) = 0). Eq. §§§ allows one 
to determine any stationary characteristic of the STO. 
For instance, the mean oscillation power p is found to be 



P 



pV (p)dp 



(10) 



Equations ([9]) and (flC)]) allow one to calculate the prob- 
ability distribution function Vo(p) and the mean mi- 
crowave power p generated by a STO in the presence 
of thermal noise and nonlinearity of the damping. The 
results of calculations performed using these equations 
are presented in Fig. [T] The main panel of this figure il- 
lustrates the probability distribution function Eq. (J5]) for 
several values of the supercriticality parameter £ = I/Ith, 
where / t h = To/ct is the critical current at which self- 
sustained oscillations start. 

In the sub-critical regime (£ < 1) the distribution 
Vq{p) has a maximum at p = 0. Expanding the in- 
tegrand in Eq. §§§ in Taylor series near p — 0, and 
keeping only the first non-zero term, one can obtain 
approximate expressions for the distribution Voip) ~ 
exp [—(1 — C)pA/(0)] and mean power in this regime, 



i-c 



k B T 



Ith - I 



(11) 



In the super-critical regime (£ > 1), the effective neg- 
ative damping T_ (p) compensates the equilibrium posi- 
tive damping T + (p) at a certain oscillator power pq. At 
this power, r_(p ) = r + (p ), and from the approximate 
Eq. ([3]) one obtains, 



Po 



C-i = I -Ith 
C + Q 1 + Q/th 



(12) 



The distribution function Vo has a maximum at p = 
Po (see Eq. ©). Thus, by expanding the integrand 
in Eq. ([9]) in a Taylor series near p — p<j, one 
can derive an approximate expression for Vo(p) ~ 
exp [— (p — po) 2 /2Ap 2 ] which is valid in the super-critical 
regime. Here Ap 2 = r ? (p )(l + Q)C/(C + Q) 2 is the level 
of power fluctuations of STO in the super-critical regime. 
The mean oscillation power p in this regime is approxi- 
mately equal to po (see Eq. (fT2"]0 . 
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Oscillator power p 

FIG. 1: Main panel: Probability distribution function Vo(p) 
(see Eq. ([9}) for several supercriticalities £. Nonlinear damp- 
ing coefficient Q = 0.3, noise level n — 0.05. Inset: Depen- 
dence of the mean oscillator power p (see Eq. $10$ ) on the 
supercriticality £ for several noise levels r/. Nonlinear damp- 
ing coefficient Q = 0.3. 

The main panel in Fig. [T] illustrates the probability 
distribution function Eq.® in the subcritical (C = 0), 
critical (£ = 1), and super-critical (£ = 2) regimes. The 
inset in Fig. [1] shows that the average power generated by 
STO (determined by Eq. (fTU)0 is substantially influenced 
by the power of thermal noise. 

As an illustration of our approach we compare the ob- 
tained theoretical results with recent measurements of 
the generated power in current-driven metallic nanopil- 
lars 0]. In the inset of Fig. [2] we show the dependence of 
the inverse mean power 1/p on the bias current 7. Ac- 
cording to Eq. (JTTJ) , this dependence is linear for small 
values of the bias current, 1/p oc (7 t h — 7). This lin- 
ear dependence in the subcritical regime 7 < 7th allows a 
method of determining precisely the threshold current 7 t h 
for microwave generation in STO for situations in which 
the influence of the thermal fluctuations is strong enough 
and determination of 2th by other means is difficult. 
Thus, from the inset of Fig.[2]one can immediately deduce 
the value of the threshold current 7th = 4.9 mA, which 
is substantially larger than the current I* « 4.2 mA, at 
which the thermally-induced oscillations become observ- 
able in the experiment 



Using the value of the threshold current 7 t h = 4.9 mA, 
derived from the experimental data in the deeply sub- 
critical regime and Eq (fT0|) . it is possible to describe the 
power generated by STO p in a wide range of bias cur- 
rents (see main panel in Fig. [2]). 

In conclusion, we have developed a simple theory of 
microwave generation in STOs based on a stochastic 
nonlinear oscillator model. We demonstrate that this 
theory allows one to accurately determine the thresh- 
old current for microwave generation from experimen- 
tal data and quantitatively describes the experimentally 
measured STO power for moderate (7 < 1.57th) values of 




Bias current / (mA) 

FIG. 2: Main panel: Dependence of the mean power p on the 
bias current I. Dots - experiment of Ref. for T = 225 K, 
solid line - theoretical dependence (|10[) for I t h = 4.9 mA, 
n = 4.2 • 10~ 4 , Q — 0.3. Inset shows the same data for inverse 
power 1/p in near-threshold range of currents. Dashed line 
corresponds to the approximate expression Eq. (|ll[l valid for 
small currents. Intersection of this line this :r-axis gives the 
value of the threshold current I t h- 



the bias current. 
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